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Temporal fluctuations in the Hadamard walk on circles are studied. A temporal standard 
deviation of probability that a quantum random walker is positive at a given site is 
introduced to manifest striking differences between quantum and classical random walks. 
An analytical expression of the temporal standard deviation on a circle with odd sites is 
shown and its asymptotic behavior is considered for large system size. In contrast with 
classical random walks, the temporal fluctuation of quantum random walks depends on 
the position and initial conditions, since temporal standard deviation of the classical 
case is zero for any site. It indicates that the temporal fluctuation of the Hadamard 
walk can be controlled. 
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1. Introduction 

Classical random walks have been used as an important tool in modern com- 
puters. It is realized by accumulating results obtained from mathematics and com- 
puter science. Quantum computer has not appeared in satisfactory form, however 
it is valuable to study both discrete and continuous quantum random walks to seek 
novel applications of quantum computers for the future. We refer to Ref. 1 and the 
references therein, concerning quantum random walks. 

In this paper we focus on discrete time case. Studies of quantum random walks 
revealed common properties and differences between classical and quantum random 
walks. If a classical random walker start from the origin on an infinite line, the 
spatial standard deviation of distribution increases in proportion to \/i after t steps. 
On the other hand, a spatial standard deviation of distribution in the Hadamard 
walk whose initial state has a value only at the origin increases in proportion to t 
after sufficiently large t steps. See Refs. 2 and 3, for examples. 

From now on we consider quantum random walks on circles. The probability 
distribution of a classical random walker on a circle with odd sites becomes uniform 
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in the limit of t — > oo, moreover time average of the probability distribution on a 
circle with both even and odd sites converges the uniform distribution. However, 
the probability distribution of quantum random walk does not converge to any 
stationary state even for odd case. On the other hand, its time average converges 
to uniform distribution for odd sites and does not converge to uniform one for even 
sites. 4 ' 5 

For that reason time-averaged distribution of quantum random walk, which 
converges a constant at a fixed site was introduced. It should be kept in our mind 
that the probability of being at position n fluctuates always bordering on time- 
averaged density. In the case of the Hadmard walk, the time-averaged distribution 
becomes uniform independently on the initial state if the number of site is odd, 
and it is agreed with the classical random walk. The fluctuation of distribution is 
not mere mathematical object, however it is related to measurement. The time- 
averaged distribution is obtained by measuring the state at a random time chosen 
in a certain interval. If the fluctuation of density function is small, we can expect 
to obtained good time-averaged distribution by the few number of measurements. 

Although details of the time-averaged distribution of the Hadamard walk on 
a circle are known, (see Refs. 4 and 5, for examples), little about the fluctuation 
of the quantum random walks is known. In this paper we introduce a temporal 
standard deviation to characterize the fluctuation of quantum random walk and 
show that the temporal standard deviation corresponding to the Hadamrd walk 
depends on the location of sites even if time-averaged density function is uniform. 
This result has a striking difference between quantum and classical cases, since 
temporal standard deviation of the classical case is zero for any site. Furthermore 
we consider the dependence of the standard deviation on the initial state. 

This paper is organized as follows. Sec. 2 gives the definition of the Hadamard 
walk and some preliminary results. In Sec. 3, we introduce a temporal standard 
deviation trjv(n) and compute it explicitly. Sec. 4 is devoted to a dependence of 
crjv(n) on the position and system size. Furthermore, Sec. 5 contains a dependence 
of (TAr(n) on initial state and parity of system size. Finally in Sec. 6 we summarize 
our results. 

2. Hadamard Walk 

Let us consider a quantum random walker who moves on cyclic sites following 
the Hadamard transform. Suppose that there are ./V sites on a circle and each site is 
labeled by an integer from and N — 1. In the Hadamard walk, the quantum state 
on a fixed site is characterized by "chirality" descried "left" and "right" . In this 
paper, the site labeled by n — 1 mod N is regarded as the site in the left to the n-th 
site. To the contrary, the site labeled by n + 1 mod TV is regarded as the site in the 
right to the n-th site. Let \L, n, t) and \R, n, t) be the wave function on n-th site at 
time t corresponding "left" and "right" , respectively. The time evolution of these 
wave functions is given by the Hadamard transform, that is, the transformation of 
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the wave function for one step is descried by 



\L,n,t) = -j=\L,n+l,t-l) + -j=\R,n + l,t-l), 

\R,n,t) = -^=\L,n- l,t- 1) - -^=\R,n- l,t- 1). 
V 2 v 2 



(1) 
(2) 



Then the total state = 0, i), |P, 0, t), |L, 1, i), |P, 1, t), ■ ■ ■ , \L, N - 

l,i),|P, TV — l,t)) T is transformed by the following unitary matrix with 2N x 
2N elements, 

P Q 

Q P 

Q P ••• 



M 



AT 



••• Q P 

Q OP 

P Q 



(3) 



where T denotes transposition and 



1 

72 



1 1 




Q = 



V2 





1 -1 



= 








(4) 



Thus the total state after t step, is given by M^^(0) for the initial state 

*(0). 

To express the total state as a function of n and t, we use known results on 
eigenvalues (Lemma 1) and eigenvectors (Lemma 2) of Mn obtained by Aharonov 
et al. 4 and Bednarska et al. 5 

Lemma 1: For j = 0, 1, . . . , N — 1 and k = 0, 1, the (2j + k + l)-th eigenvalue 
corresponding to Mm is given by 



Cjk — 



V2 



(-1)* 



cos^ 



2J7T 

N 



% sin 



2J7T 



We here define v° k [ and w| fcjZ to express eigenvectors by 

v j.k,i = ajkbjk^N, 
v lk,i = ajk^N, 



where 



ojn = e N 
a 3 k = 



v(i + IVI 2 ) 



b jk = u ] N { (-If a/1 + cos 2 + cos£,- 



(5) 

(6) 
(7) 

(8) 
(9) 

(10) 
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= 



2j7T 



(11) 



* N 

Lemma 2: For Z = 1,2,..., 2N, the l-th element of the eigenvectors corresponding 



to Cjk is given by 



v 3,k,l 



V lk,(l+l)/2 ( l = ° dd )' 



(12) 



I V 3-,k,l/2 



(I = even), 

For the convenience of readers, from now on here we check directly that eigen- 
vectors Vj t k — ( v j.k,i, ■ ■ ■ ,Vj,k,2N) T satisfy the next equation 

M N v jik = c jt kVj,k- (13) 

(a) Suppose I is odd, then the (2m — l)-th element of the left-hand side of the 
Eq. (13) for m = 1, 2, • • • , N is given by 

(2m - l)-th element of LHS 
— ~~J^ (^j.fe.m+i mod N + v j,k,m+i mod n) ' 

2j(m+l)i 



V2 
djk 
V2 



1 + e— (-1) \/ 1 + cos 2 & + cos& 



(14) 



where ^ = 2nj/N. On the other hand, the right-hand side of Eq. (13) is given by 
(2m - l)-th element of RHS 

= CjfcajkbjkLO^ , 

CLjk 



V2 



2j(m+l)j 



(-l) fc ,/l + cos 2 ^ + cos£. 



x ( (-l) fe y / l + cos 2 ^- +isin£j 



(15) 



(b) Similarly suppose I is even, then the 2m-th element of the left-hand side of 
the Eq. (13) for m = 1, 2, • • • , N is given by 

2m-th element of LHS 

= — (v° -v e \ 

fn \ j,k,m-imodN j,k,m—l mod N J ' 



V2 

cijk 

71 



2 3 '(m-l)«7r 




2«H / 


e « 


-e « 1 









(-l)V 1 + cos 2 ^- + cos£j 



(16) 



On the other hand, the right-hand side of the Eq. (13) is given by 

2m-th element of RHS 

_ j m 

— Cjkdjk^N ) 



djk 

V2 



e 



(-l) fe y / l _ fco^ 2 C~ + isin£. 



(17) 



Using the Eular's formula, we can check that both sides of Eq.(13) coincide for each 
case. 
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3. Temporal Standard Deviation 

Since we obtained the eigenvalues and eigenvectors, matrix Mm is transformed 
into a diagonal matrix and the wave function after t steps is generally expressed by 
a liner combination of c* fc . In real system we can not observe the wave function but 
probability of being at position n at time t. This probability P/v(^, t) is calculated 
from the wave functions by 

P N {n, t) = (L, n, t\L, n, t) + (R, n, t\R, n, t). (18) 

In contrast to classical random walks, the probability P/v(n, t) dose not converge 
in the limit of t — > oo. Thus we consider a time- averaged distribution defined by 

1 T_1 

P N {n) = lim - VPiv(rM), (19) 

if the right-hand side of Eq. (19) exists. Aharonov et al. 4 showed that Pv(n) exists, 
and it is independent of both an initial state and n if all eigenvalues of Mm are 
distinct. If the number of site is odd, then all eigenvalues of Mm are distinct, so we 
immediately have Pjv(ft) = 1/N for any n = 0, 1, . . . ,N — 1. On the other hand, 
if N is even, then there exist degenerate eigenvalues of Mm and it is possible to 
derive non- uniform distributions. In the case of the classical random walk, the time- 
averaged distribution of a walker becomes uniform independently on the parity of 
the system size. Thus we can not distinguish the classical random walk and the 
Hadamard walk on a circle with odd sites only by the time-averaged distribution. 

As we mentioned above, the P/v(n, f) does not converge in the limiting t — > oo. 
It means that the probability always fluctuates to the upper and lower sides of 
Pjv(n)- For this reason, we define the following temporal standard deviation o m(ji)'- 



a N (n) = 



\ 



1 T_1 

lim -]T(Pv(M)-PvW) 2 , (20) 

1 — »oo 1 — 



t=0 



if the right-hand side of Eq. (20) exists. 

Here we consider the classical case. In the case of classical a random walk starting 
from a site for odd sites (i.e., aperiodic case), there exist a € (0, 1) and C > (are 
independent of n and t) such that 

\P N {n,t)-P N (n)\<Ca\ 

where Pv(^) = 1/A for any n = 0, 1, . . . , N — 1, (see page 63 of Schinazi, 6 for 
example). Therefore we obtain 

IX:(P w (n,t)-P w (n)) 2 <^i-J. 

t=o 

The above inequality implies that for any n = 0, 1, . . . , N — 1, 

a N (n) = 0. 
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in the classical case. As for N — even (i.e., periodic) case, we have the same con- 
clusion (7jv(n) = for any n by using a little modified argument. 

In this situation, we first ask a natural question whether the ctjv(ti) is always 
zero or not in quantum case. Moreover if a\/v(n) is not always zero, then we next 
ask a question whether the <JN{n) is uniform or not. 

From now on we concentrate our attention to the system with odd sites. Fur- 
thermore we assume that the initial state is a fixed ^(0) = (1,0,0,- • - ,0) T . Let 
us express the standard deviation CTjv(n) starting *(0) = (1,0,0, • • • ,0) T by wave 
functions. To do so, plugging Eq. (18) into Eq.(20) we have 

a%(n)=lim^-Y,((L,n,t\L,n,t) + (R,n,t\R,n,t}--j , (21) 
t=o ^ ' 

where we used the result Pin) = 1/N for odd N. Since the matrix Mm was already 
diagonalzcd in the previous section, the wave function \L,n,t) and \R,n,t) arc 
easily obtained as 

N-l 

\L, n,t)=Y^ (^oc*o + aj-icJ-O , (22) 

3=0 
N-l 

\R,n,t) = (Mo+Mi)- (23) 

3=0 



where 



(-l) fc e 



v /l + cos 2 ( ^)^ 



&k = v (25) 

2^1+0082(2^) 

Since the matrix Mjv is unitary matrix, the eigenvalue Cjk is written as e l ° ik where 
6jk is the argument of Cjk- Thus the probability (L,n,t\L,n,t) is expressed by 

JV-l 

(L,n,t\L,n,t) = ^ (|a, | 2 + Ki| 2 ) 

3=0 

N-l 1 

+ E E ^•i,fcofci«jofco«; i fc 1 e <( ^*o-«ii*i)* ) (26) 

30,31=0 fc o ,fci=0 

where 

c _ / jo = ji and k = ki , . 

d ^> fe ° fe i-\l otherwise. 10 

The first term of Eq. (26) is a constant and the second term vanishes by ap- 
plying an operation liniT^oo yX^tlo 1 - Similarly we can divide (R, n, t\R, n, t) into 
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a constant term and a vanishing term: 

JV-l 

(R,n,t\R,n,t) = ]T (|/? j0 | 2 + \Pji\ 2 ) 



3=0 

JV-l 1 

+ E E ^.fcofc^iofco^fe^o*"-^* 15 *. (28) 
jo ,ji=0 fe ,fei=0 



As a result we have 



JV-l 

P(n) = (l«,o| 2 + M 2 + Ift-ol 2 + l/3,i| 2 )- (29) 

Plugging Eqs. (24)-(25) into Eq. (29) gives confirm P(n) = 1/JV for any n — 
0,1,..., JV-l. 

Let express <J 2 N (n) as function of eigenvalues. By using Eqs. (21), (26) and (28), 
we get 

^ T-l JV-l 1 

t=0 h,ji,h,33=0 ka,kite,k 3 =0 
x(«iofco«i lfcl + ft»i,^ ltl )(%fe4i3 + ^fei^fe)^, (30) 

where 

A0 = Ad(j ,k ,j 1 ,k 1 ,j 2 ,k2,j3,h) 

= Qjoko - 8jiki + 6j2k 2 - 9j 3 k 3 - ( 31 ) 

If A6> 7^ (mod 27r), then hmr^oo e lAet /T converges to a zero. 

Thus the variance Ojy(n) is obtained by taking summation of (ctjoko^ifci + 
/ 3 jokoPj 1 k 1 )( a j2k 2 a* 3k3 + P j2 k 2 P* 3k3 ) over combinations j , k , ■ ■ • , j 3 , fc 3 satisfying 
A# = (mod 27r). For this reason, we consider the combinations which satisfy 
AO = (mod 2n) for a given combination (jo,ko,j\,k\). We define the following 
four conditions: 

(a) Re[cj fe ] = Re[cj lfel ] and Im[c Jofeo ] = -lm.[c jlkl ], 

(b) Re[cj fe ] = -Re[c jlkl ] and lm[c joko ] = lm[c jlkl ], 

(c) Re[c joko ] = -Rc[c jlfel ] and Im[c jofeo ] = -Im[cj lfcl ], 

(d) Rc[c joko ] = Re[c jlfel ] and Im[c Jofeo ] = Im[c Jlfel ]. (32) 

If each condition from (a) to (d) is not satisfied, there are four different combinations 
which satisfy A9 = (mod 2tt) for a given combination (jo, ko,ji,ki): 

(a) j2=N-j mod JV, k 2 = k , j 3 = N — j 1 mod JV, fc 3 = fci, 

(b) J2 = jo, fc 2 = 1 - fc , j 3 = ji, k 3 = l- fci, 

(c) .72 = JV - ji mod JV, fc 2 = 1 - fci, .73 = JV - jo mod JV, fc 3 = 1 - fc , 

(d) j 2 = ji, k 2 = h, j 3 = j , &3 = feo- (33) 
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If one of condition in (32) is satisfied, we find the same values in (33). 

We can introduce a geometrical picture into these combinations. All eigenvalues 
exist on a unit circle in a complex plane and there is a symmetry in these points. For 
example, when we find Cj k at the point A, we can find always another eigenvalue 
on a point which is obtained by reflecting point A in the x— axis and its evienvalues 
is expressed by c n _. q moc [ N ko . We also find always an eigenvalue on a point which 
is obtained by reflecting point A in the y— axis and its eigenvalues is expressed 
by Cj 0; i_fc . Additionally we note that if the relation between Cj a k a and Cj 1 k 1 is 
symmetrical with respect to the origin and the relation between Cj 2 k 2 and Cj 3 k 3 is 
also symmetrical with respect to the origin, then A9 = (mod 2w) . 

From the relation (33) we express J2,fc2j3 and k% as a function of jo, fco, ji 
and k\. Thus cr N {n) is obtained by summations over jo, ko, j\ and k\. Taking the 
degenerate into account correctly we obtain a%(n) by some computations. Before 
showing the result we define next functions 

JV-l 

*o = E jt^t, (34) 

3=0 3 



3=0 
N-l 



^ cos((ra-l)flj) + cos (nfy) 
S+{n > ^ 3 + costf, ' (36) 

3=0 J 

c _ v- 1 cos (( n z ])M z cos ( ng j) cif\ 



3 + cos Oj 

3=0 J 
N-l 



" 1 7 + co S 20 J -+8co S J - C os 2 [(n-l)0 J -] 

S2in) = IrTT^TVi ' ( 38 ) 



(3 + cos^) 2 

where 9j = 4nj/N. Now we show cr^(n) as a main result in this paper: 



Theorem 3: When A is odd, we have 



= ^4 t 2 + + US o + 10S oSi + 3Sf - S 2 (n)] - A (39) 

for any n = 0, 1, . . . , N — 1. 



4. Dependence of crisr{ri) on the position and system size 

The obtained formula for ajsi{n) is written in the form of a single summation, 
therefore, it is calculated very easily in comparison with Eq. (30) containing nine- 
fold summations. We show <7jv(n) for A = 3, 5, • • • , II in Fig. 1. 

First we find that crjv(n) depends on the position and it has the maximum value 
at n = and 1. The dependence of <rjy(n) on the position from 2 to A — 1 is not 
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0.3 
0.25 
^ 0.2 
b 0.15 
0.1 
0.05 

0123456789 10 
n 

Fig. 1. Dependence of cr^(n) on the position and system size 

clear only from Fig. 1, however we find that cr/v(n)= (Jn(N + 1 — n mod N) for any 
n and all values are distinct except a pair n and (N + 1 — n) mod N. 

Next we consider asymptotic behavior of ctjv(0) for large N to observe the 
dependence on the system size. Setting n = in Eq.(39) gives 




El cos Oj 
J - 2 
3 + cos 9 j 

3=0 J 




' ^ 15-11 cos 6j\ ^9 + 4 cos 6j + 3 cos 26 3 
3 + cos0j I ^ (3 + cos^) 2 



On the other hand, we have 

1^1 I / iT I I 

~r~ I t. dx = — =, 

47T Jo 3 + cosx 2V2 



1 jV_1 1 

um - y — 

n^oo N 4^ 3 + cos 6 

J=0 



1 ^-i A 

1 v-^ COS ft, 
hm — > — = 

3=0 



r 47r i 
:/ 3T^ fl 

cosx , 3 
ax = 1 = , 

> 3 + cosx 2V2 



lim 



1 JV " 1 1 

..... 1 V - 

jv^oo N 4-i (3 + cos 9 7 -) 2 

1 

lim — >^ 

AT^oo N ^ (3 + COS 0,- 



47T J a 



>4tt 



COS 0j 



1 



/■47T 



) 2 47T 7 

1 



•47T 



lim 1 V C0S(2 ^ } - 

v^oo N (3 + cos^) 2 4tt7 (3 + cosx) 2 



(3 


+ cosx) 2 




cosx 


(3 


+ cosx) 2 




cos 2x 



1 

16V2' 
45 



(40) 

(41) 
(42) 

(43) 
(44) 
(45) 
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Using the above results, we obtain asymptotic behavior of the variance cr^(O) for 
sufficiently large N as follows: 



Proposition 4: For N — > oo, we have 



13-8V2 7\/2-16 



N 2 



■ + 



2N 3 



+ 



1 

N3 



(46) 



The above result implies that the fluctuation cxjv(0) decays in the form 1/N as 
N increases. Fig. 2 shows the comparison between exact values of ctjv(O) and the 
approximate values obtained from Eq. (46). We see good agreement between them 
except N = 3. 



0.3 
0.25 



o 0.2 



z 
b 



0.15 
0.1 



10 

N 



15 



20 



Fig. 2. Comparison exact values of <Tjv(0) (solid circles) with approximate values (solid line) 
given by Eq.(46) for N = 3, 5, • ■ • , 21. 



5. Dependence of crjv(n) on Initial States and Parity of System 
Size 

5.1. Dependence of cr N (n) on Initial States 

In the previous sections, we have shown the results covering only at the case 
where the system size is odd and the initial state is \L, 0, 0) = 1. The time averaged 
distribution on a circle including odd sites is independent of the initial state. Is 
also the temporal standard deviation independent of the initial state? Since the 
temporal standard deviation depends on the position, it is clear that the temporal 
standard deviation depends on the initial state. We consider here whether the tem- 
poral standard deviation depends on the initial state even if the initial probability 
distribution is the same. 

In order to observe the dependence of the temporal standard deviation on initial 
state without changing the initial probability distribution, we set |L, 0,0)=a and 
\R, 0, 0)=-^! — a 2 i where a € [0,1]. As a particular case, we have a symmetrical 
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initial state by setting a = l/y/2. Thus we can observe the dependence of ffN(n, a) 
on symmetry break down by changing a. We first present the result for N = 3. The 
standard deviation <r 3 (n,a) for the above initial state is obtained as follows: 

<73(0,a) = ^p, (47) 

a 3 (l,a) = -^ v / 96o! 4 -75a 2 + 25, (48) 
45 

<7 3 (2,a) = — v/96a 4 - 117a 2 +46. (49) 
45 

Fig. 3(a) shows the dependence of 03(71, a) on the parameter a. One clearly finds 
that (73(0,0;) is independent of the initial and both <r 3 (l,a) and 173(2, a) depend 
on a. We further find that 03(1,0;) is equal to o 3 (2,a) at the symmetrical case 
a = 1/V% and the minimum values are located near a — l/y/2. 

Fig. 3(b) and Fig. 3(c) show approximate values by setting T = 10 4 in Eq. (19) 
and Eq. (20) in place of taking a limit for N = 5 and N = 7. Similarly with Fig. 
3(a), the temporal standard deviations at the position n = are independent of the 
parameter a and the minimum values arc located near a = 1/ V2. The dependence 
of ojv(n, a) on the position is weak between 2 and N — 2, and they seems to be in 
close on a single line. 

5.2. Temporal Standard Deviation on a Circle with Even Sites 

We have concentrated our attention into the case where the number of the system 
size is odd. The main reason why we avoid even cases is that the eigenvalues of 
matrix Mjv are highly degenerate, and the time averaged distribution itself depends 
on the position of site and the initial state. While the eigenvalues for even cases 
is exactly obtained by Eq.(5), there are many possible combinations of eigenvalue 
which satisfy A6 = (mod 2n) for even case. Thus we can not calculate ojv(n) 
in the same way with odd cases. Therefore, we here try to carry out approximate 
calculations instead of seeking analytic results. Fig. 4 shows the approximate ojv(n) 
for N = 4, 6, • • • , 12 under the same initial condition \L, 0, 0) = 1 and T = 10 4 . 

Let us compare Fig. 1 with Fig. 4. First we find that the values of the temporal 
standard deviation for fixed N is almost the same except n = 0, while the values 
for even N is strongly dependent on the location of the site in comparison with 
odd cases. Second we find that the maximum values of ojv(«) for odd N exist 
only at n — and n = 1 as shown Fig. 1, while the four same peaks are found for 
N = 4,8, 12 in Fig. 4. We confirm numerically that this behavior is always observed 
up to N — 22 when N is a multiple of 4. 

6. Summary 

We expressed the temporal standard deviation of the Hadamard walk on a cir- 
cle with odd sites for a pure initial state \L, 0, 0) = 1 in an analytical form. The 



12 




formula is obtained by calculating the production between four Fourier coefficients 
which come from production of "L" and "i?" states. Unlike the time-averaged dis- 
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Fig. 4. Dependence of &pf(n) on the position with even system size. 



tribution, the temporal standard deviation depends on the location of the site. The 
fluctuations decrease in inverse proportion to the system size for large N and its 
coefficient was calculated rigorously. 

When we set the initial state in the form \L, 0, 0)=a and \R, 0, 0}=\/i — a 2 i, the 
analytical result for N = 3 and numerical simulations indicated that the temporal 
standard deviation depends on the initial state except n — 0. We would speculate 
that the temporal standard deviation takes the maximum value at n — and its 
value is independent of a. The values of the temporal standard deviation at neighbor 
site of origin is somewhat larger than other sites and the remain are almost the 
same. The dependence of the temporal standard deviation on the position for even 
system size is much complex in comparing with that for odd system size due to 
degeneration of eigenvalues. Several peaks are observed when the system size is a 
multiple of 4. 

From the viewpoint of technology, it seen to be of value to consider whether the 
temporal standard deviation is controlled or not. It is shown that the time-averaged 
distribution can not be controlled for odd system size, however, for even system size, 
it is controlled thanks to the degenerate of eigenvalues. Our results show that the 
standard deviation can be controlled. First the dependence of temporal standard 
deviation on location is used. Second the dependence on initial state is also used. 
However we note here that ctjv(O) for N = 3 is independent of the initial state. In 
this meaning, its control is limited. 
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